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Abstract 

We consider Angelesco ensembles with respect to two modified 
Jacobi weights on touching intervals [a,0] and [0,1], for a < 0. As 
a —7- —1 the particles around experience a phase transition. This 
I '"'i transition is studied in a double scaling limit, where we let the num- 

ber of particles of the ensemble tend to infinity while the parameter a 
1^ tends to —1 at a rate of 0{n~^''^). The correlation kernel converges, 

Q> in this regime, to a new kind of universal kernel, the Angelesco kernel 

■^ K ^^. The result follows from the Deift/Zhou steepest descent anal- 

og ysis, applied to the Riemann-Hilbert problem for multiple orthogonal 

^J polynomials. 

o 

Tl 1 Introduction and statement of results 

K> Multiple orthogonal polynomial (MOP) ensembles [21j form an extension of 

H the more familiar orthogonal polynomial (OP) ensembles [20]. The latter 

appear as the eigenvalue distributions of unitary random matrix ensembles. 
A major topic of interest in random matrix theory deals with the uni- 
versality conjecture. This conjecture claims that, as the size of the matrices 
tends to infinity, the local eigenvalue statistics of a random matrix ensemble 
converge to universal limits, i.e., limits that are independent of the precise 
probability distribution on the matrices. The limits only depend on macro- 
scopic eigenvalue characteristics and on the symmetry class of the matrix 
ensemble. 



For the class of unitary random matrix ensembles the problem concern- 
ing universality can be translated into a problem concerning the associated 
sequence of orthogonal polynomials (OP). By the Gaudin-Mehta method 
[TTl [HI, [25] the eigenvalue density function can be written as a determinantal 
point process. The relevant correlation kernel is the OP kernel, which is also 
known as the Christoffel-Darboux kernel. Universality is then expressed as 
the convergence of this OP kernel to a certain universal limiting kernel. This 
convergence can be derived by considering the large degree asymptotics of 
the OP. Common limiting kernels include the sine, Airy and Bessel kernels. 

For other random matrix ensembles a connection exists with multiple 
orthogonal polynomials (MOP): the eigenvalues are a determinantal point 
process with the so-called MOP kernel as correlation kernel. A prime example 
of such a connection occurs in the unitary random matrix model with external 
source. By analysis of the MOP kernel in this model universality has been 
shown and the sine and Airy kernels appeared [2l |5] . However, in a critical 
case this model also exhibits behavior that does not appear in OP ensembles: 
under certain scaling the MOP kernel converges to the so-called Pearcey 
kernel [6l [HI [29] . MOP ensembles also gave rise to new critical behavior in 
the papers [Til 123]. 

This is a recurring feature: MOP ensembles can exhibit a wider variety 
of critical phenomena than OP ensembles. In this paper we consider another 
type of MOP ensembles, the so-called Angelesco ensembles. These ensembles 
do not appear in a natural way as the eigenvalue process of a random matrix 
ensemble, but the notion of universality still makes sense. We show that a 
new kind of critical behavior occurs in Angelesco systems, determined by a 
new kind of universal kernel, that we call the Angelesco kernel K^'^^. 

An Angelesco system of weights P, [2S] is a system of r > 2 weights 
w = {wi, . . . ,Wr) on the real line such that the supports of the weights 
are contained in intervals with pairwise disjoint interiors. Thus there exist 
bounded real intervals Ai, . . . , A^ such that 

AjnAk=0 ioij^k (1.1) 

and Wj is non-negative and integrable on the real line with 

Wj{x)=0 ioTxeR\Aj. (1.2) 

Angelesco ensembles are MOP ensembles associated with Angelesco sys- 
tems of weights. Let n be a multi-index (ni, . . . ,nr) G N'' and define |n| := 



Hi + ■ ■ ■ +nr- The Angelesco ensemble with respect to w and n is the prob- 
abihty measure on RI"' with probabihty density function 



p^(xi, . . . , xiH\) ■■= y: d^* K ] i, fc=i ■ ^^* [/i(^fc)]rfe=i 



i|n| 



where Z^ is a normahzation constant and 



Jn\~ 



-"i-i+J 



X 



x-' ^Wi(x) 



for j = l,...,ni,z = l,...,r. 



:i.3) 



:i.4) 



The Angelesco properties ( 1. 1 )-( 1.2 ) imply that ( 1.3 ) is non-negative for every 



la;i, 



X 



\n\, 



G M'"' , and can be non-zero only if n,- of the particles X\ 



,x 



\n\ 



are in Aj for every j = 1, . . . , r. Thus in an Angelesco ensemble, with 
probability 1, n,- particles are located in A,- for every j = 1, . . . , r. 



:i.5) 



The probability density function (1.3) is a biorthogonal ensemble [7], 
which is a special case of a determinantal point process. The correlation 
kernel is 

\n\ \n\ 

j=i fc=i 
where {M~^)kj is the kj-th entry of the inverse of the matrix 

M = {mj^k)^^l^i , rrij^k = / x^~^ fk{x)dx. 



:i.6) 



|n| 



\n\\ 



dei[Kn{xj,XkWL^ 



That is, see m EI]: 

and for each m = 1, . . . , |rr|. 



:i.7) 



:i.8) 



where pm is the m-point correlation function. Another representation for Kf^ 
is by means of a determinant 



Knix,y) 



detM 



M 



j.|H|-l 

/i(y) /2(y) - f\n\(y) 
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Figure 1: The phase diagram for Angelesco ensembles with two weights, 
supported on [a, b] and [0, 1]. The gray areas represent values of a and b that 
are forbidden because of the restriction a < b < 0. Remark that the scale on 
the two axes is different. 



where M is the moment matrix (1.6). 



In this paper we consider Angelesco ensembles with respect to two weights 
Wi and W2- Also we restrict ourselves to diagonal multi-indices n = {n,n). 
Our interest is in the local asymptotics of the MOP kernel Ka = Kn,n as- 
sociated to the weights wi and W2- By rescaling if necessary we can assume 
that Ai = [a, b] and A2 = [0, 1] for some a < 5 < 0. 

The type of local behavior around a certain point is suggested by the 
behavior of the limiting mean particle density in that point. Assume that 
Wj > a.e. on Aj. Then this limiting density only depends on the endpoints 
of the two intervals Ai and A2. Depending on the values of a and b we get 
a qualitatively different picture. We can distinguish 3 cases and a number of 
phase transitions, see Figure [T] 

Case III corresponds to relatively large separation between the two inter- 
vals. As a consequence there are hard edges in all four endpoints a, b, and 
1, meaning that the limiting particle distribution blows up like an inverse 
square root. After rescaling around one of these endpoints the MOP kernels 
converge to limiting kernels that depend on the behavior of the weight at 
that endpoint. In the typical case of Jacob i-type weight functions 



wi{x) = (x - a)"(6 - x)'^hi{x)x[a,b]i^): 
W2{x) = x^{l - xYh2{x)x[o,i]ix), 



1.9) 



this leads to Bessel kernels. 

In the cases I and II the separation is small. A gap will emerge on the 
larger of the two intervals Ai and A2, due to the pushing of the zeros effect. 
There will then be three hard edges, and one soft edge, where the limiting 
particle density vanishes like a square root. Around the soft edge the rescaled 



MOP kernel will converge to the Airy kernel (in the case of weights (1.9)). 

The transition between case III and cases I and II is probably related to 
the Painleve II equation as in [9]. The Bessel, Airy and Painleve II kernels 
already appear as scaling limits for OP ensembles [T5] . 

This paper deals with the critical point (a, 6) = (—1,0), where the two 
intervals Ai and A2 are touching and of equal size. In this case the limiting 
particle density behaves like |x|~3 as x — )■ 0, see [13], [19]. We study this 
critical case in a double scaling limit, where we put 6 = and let a tend to 
— 1 as n — )■ 00: 

\/2t 
a = a„:=-l + ^^ for r G M. (1.10) 

The parameter r is known as the double scaling parameter. This scaling 
regime corresponds to the phase transition between cases I and II with b 
fixed at 0. 

For a 7^ —1 and 6 = we have a gap in the spectrum between and s 
where s = s^ is given by 

^ - ''' + ^'' (1,11) 



" 9(a2-a + l)' 

see [19]. In s a soft edge appears, while in we have a hard edge. Plots of the 
limiting particle densities in the cases a<— 1, a = — 1 and a > —1 are given 
in Figure |2j As a tends to —1 the gap closes and the soft edge meets with 
the hard edge. This soft-to-hard-edge collision will give rise to new critical 
behavior in 0. 

So for the remainder of the paper we take for a < 

Ai:=[a,0], A2:=[0,l], (1.12) 

and as in [13] we consider the following modified Jacobi weights: 

wi{x) := {x - a)''\xfhi{x)x[a,o]{x), 
W2{x) := a;^(l - xyh2{x)x[o,i]ix), 



rL 





so 

a < -1 

(a = -1.5) 





a = —1 




Figure 2: Plot of the limiting particle densities in Angelesco ensembles with 
two weights on [a, 0] and [0, 1] in the cases a < —1, a = —1, and — 1 < a < 0. 
For clarity the size of the gap between and s has been exaggerated in the 
figure. In the left and right plots the density has a local maximum near 
s. Even though the plots may suggest otherwise, the density is in fact real 
analytic there. 



with q;,/3,7 > —1 and functions hj that are positive on A^ and analytic in 
a neighborhood of A for j = 1,2. Note that we use the same exponent /3 in 
both wi and W2. 



In [13] we considered the weights (1.13) and obtained a Mehler- Heine 
asymptotic formula for the MOP associated to wi and W2 and n = (n, n) as 
n^ oo. See also [27] 128]. 

We use the notation Kn^n{-, S ^t) to denote the dependence of the correla- 



tion kernel (1.5) on the parameter a. The main result is then as follows: 



Theorem 1.1. Let Kn,n be the correlation kernel of the MOP ensemble with 
weights (1.13) and diagonal multi-index n = {n,n). Let T,x,y G M with 



y 7^ 0. Then as n ^ +oo we have: 



V2: 



—K 

3 ^^n,n 



X 



na 



V2nr 



V2'. 



-;-i + 



n2 



V2i 
na 



K^'^^{x, y;T) + 



1 



(1.14) 



for some limiting kernel K^°s(x, y;r) that will be described below in Propo- 
sitions 1.2 and 1.3. The O-term in (1.14) is uniform for x in a bounded set 



and y in a compact subset m M \ {0}. 

The limiting kernel K^'^^ is referred to as the Angelesco kernel. It has an 
expression in terms of the Angelesco model parametrix \1' which we will de- 
scribe in Section [2^ More explicit expressions exist, one involving a certain 
pairing of two analytic functions qo, tq and one in terms of a double contour 
integral. 
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Figure 3: The contours Fq and Fq appearing in the integral formulas (1.15), 



(1.16), and (1.20) for the Angelesco kernel K^'^^. The dashed line denotes 



the branch cut of t^ and s'^ 



Proposition 1.2. Define functions go one? tq by 

1 



qo{x) = qo{x;T) := ^^-\x 



ro{y) =ro{y;T) := —^\y 
Zm 



27rz jp| 

1 



t~^"3e~~i^+*rft 



-/3-1 



s^e- 



xeR, (1.15) 

ds, yeR\{0}. (1.16) 



ilfere the contours Fq one? Fq are shown in Figure\^ and we choose the prin- 
cipal branches for the fractional powers 1 1— )■ t~^~^ and s ^ s^ . Then we can 
write K^'^s(x, y; r) as a pairing of qo{x) and ro^y) in the following way: 



K^^^ix, y;T) 



yroiy)Qoix) - ((/3 + ^)ro{y) + yr'oiy)) q'oix) 

+ iyr'oiy) + W + 2)r[,(2/) - Tro{y)) qo{x) 



(1.17) 



From the expressions for go and tq it is clear that we have 

go(x;r) = go(-x;-r), ro(i/; r) = ro(-|/; -r). (1.18) 

Hence the Angelesco kernel satisfies the following (expected) symmetry: 



"Ane 



{x,y;r)=K^-^{-x,-y;-T). (1.19) 

Proposition 1.3. Let the contours Tq and Fq be again as in Figure\^ and 
take x,y G M \ {0}, r G M. Then we have the following double integral 
formula for K'^'^s(x, y; r) .■ 



K^°^(x,y;r) 



sgn(y) 



TX __ X 

e * 2iP 



+t 



ds dt. 



1.20) 



{27iiY Jt^rjs^v,t^xs-yt^^-l 
Again main branches are used for the fractional powers t ^-^ t^ and s \-^ s^ . 




Figure 4: Numerical simulation of 50 non-intersecting squared Bessel paths. 
On the horizontal axis, time is running from to 1. The vertical axis has 
the position x. The figure is taken from 



Remark 1.4. In [22], [23] the authors consider a model of n non- intersecting 
squared Bessel paths, conditioned to start at time t = in the point x = 
a > 0, and to end at time t = 1 at x = 0, see Figure |4j For each fixed time 
t the particles form a MOP ensemble with respect to orthogonality weights 
involving the modified Bessel functions J/j and //3+1 for some (3 > —1. 

Associated to each fixed time t one has a limiting mean particle distri- 
bution as n — )■ 00. For small values of t the particles are supported on an 
interval in M+, away from 0. On the other hand, for values of t near 1 the 
point is contained in the support, meaning that particles stick to the hard 
wall in 0. There is then a critical time t = t* & (0, 1) at which the gap be- 
tween the particles and the wall in closes. This phase transition is similar 
to the phase transition in Angelesco ensembles. 

By taking a double scaling limit around a; = and t = t* a new kind 
of universal limiting kernel, K'^"^, was obtained ^B] equation (1.18)]. The 
function K'^"^ depends on two positive position variables x and y, a double 
scaling parameter r, and the parameter /3 appearing in the modified Bessel 
functions. It turns out that this critical kernel is almost the same as the An- 
gelesco kernel K^°^. By applying a simple substitution to the double integral 



formula for K'^^^^ [231 equation (1.19)] one obtains the following identity: 

K^''\x,y;T) = ^Kf\y,x-T), x,y > 0, (1.21) 

where we used the notation Ko"^^ to indicate the dependence of the kernel 
on the parameter /3. 

Although the kernels are different this identity implies an equality at the 
level of the correlation functions. Hence, after proper identification of the 
parameters, the critical behavior in Angelesco ensembles coincides with the 
critical behavior in the model of non-intersecting squared Bessel paths. A 
significant difference between the two models is that in the squared Bessel 
paths model the particles remain on the positive half-line, while in our model 
they are located on both sides of and there is interaction between the two 
groups. 



The proof of Theorem 1.1 follows from the application of the Deift/Zhou 



steepest descent method on the Riemann-Hilbert (RH) problem for MOPs. 



In Section^ we recall the RH formulation for the weights Wi and W2 (1.13) 
and we give a formula for the MOP kernel in terms of the RH problem 
for MOPs. Also we recall the RH problem for the Angelesco local model 
parametrix \1/ which is taken from [13] . 

The Deift/Zhou method consists of a sequence of invertible transforma- 
tions reducing the original RH problem into a normalized RH problem, for 
which uniform estimates can be made. The analysis is the same as in jT3] . 
Instead of giving full details we will, in Section |3| give a quick overview of 
the transformations and describe their effect on the expression for the MOP 
kernel. 

In Section |4] we shall use the series of transformations to show that the 



limit in Theorem 1.1 indeed holds, and we give expressions for the Angelesco 
kernel ]K^°^ in terms of the Angelesco parametrix \l/. 

Finally, Section [5] contains the proofs of the two explicit expressions for 



K^'^s stated in Proposition 1.2 and Proposition 1.3 



2 Riemann-Hilbert problem 

This section recalls the Riemann-Hilbert (RH) characterization of the rele- 
vant multiple orthogonal polynomials and the Angelesco local model parametrix. 



2.1 Multiple orthogonal polynomials 

We state here the RH problem for MOPs [30] with respect to the system 



of weights {wi,W2) given in (1.13) and the multi-index {ni,n2) € N . The 



endpoint conditions are as in ^^. 

We look for a 3 x 3 matrix- valued function Y such that: 

• y is defined and analytic on C \ [a, 1] . 

• Y has continuous boundary values Y± on (a, 0) and (0, 1) and they 
satisfy the jump relation Y^ = Y^Jy with Jy given by 

(1 Wi{x) W2{x)\ 
1 \ , xe (a,0)U(0,l). (2.1) 

1 ; 

Recall that wi{x) = outside of [a, 0] and ^2(3;) = outside of [0, 1]. 

• As 2; —> 00 we have 

\ni-i-n2 g g 

Y{z)= [I + 0[-] ] \ ^-"1 I . (2.2) 





Z-' 

Y has the following behavior at the endpoints of the intervals: 

/I e{z) 1\ 
Y(z) = \1 e(z) 1 1 , as z ^ a, 

'{z-a)" ifa<0, 

where e{z) = ^ \og{z — a) if a = 0, (2.3) 

1 if a > 0. 



Yiz) 




[z-1)^ if7<0, 
where e{z) = { \og{z - 1) if 7 = 0, (2.4) 

if 7 > 0. 
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Y(z) = \ 1 e(z) e(z) \ , as z ^ 0, 




if /3 < 0, 
where e{z) = ■{logz if /3 = 0, (2.5) 

1 if /3 > 0. 

The O-symbol is to be taken entry-wise. 

This RH problem has a unique solution in terms of MOP with respect to 
the weights wi and W2. In particular we have that Yii{z) is a monic polyno- 
mial of degree ni + n2 that satisfies the multiple orthogonality conditions 

x^Yu{x)wj{x)dx = 0, k = 0, . . . ,nj — 1, j = 1,2. 

For the full expression of Y we refer to [13] . Most important for the present 
paper is the following formula for the MOP correlation kernel in terms of Y, 
see Ullin]: 




^'^^•"^^'^'^^ = 2m{x-y) ^° "^'^^^ "^'^^^^ F+(z/)-^F+(x) | | . (2.6) 

This formula is valid for x 7^ y G M with y ^ {a, 0,1}. As before it is 
understood that Wi{y) = for y outside (a, 0), and 1^2(2/) = ioi y outside 



(0, 1). The formula (2.6) can be extended to the case x = y ^ {a, 0, 1} by 
I'Hopital's rule. 

2.2 Angelesco local model parametrix 

At a crucial step in the steepest descent analysis we need to do a local 
analysis at the point 0. This step involves certain special functions, which 
are combined into a 3 x 3 matrix-valued function, called the Angelesco model 
parametrix \l/ and which was introduced in fiSl Section 2.2]. The function 
\l/ depends on two parameters (3 > —1 and r e M. Since (3 is considered 
fixed we do not emphasize the dependence on /3. We may write \&(2;;r) to 
emphasize the dependence on r. 

The function \1/ has the following RH characterization: 

11 



aP-^i 




Figure 5: The contour E,j, and the jump matrices of \l/. 



^ is defined and analytic on C \ S^ where S^ is a contour consisting 
of six oriented rays through the origin, as shown in Figure [5j 

\1/ has continuous boundary values on S,j, \ {0} that satisfy the jump 
condition 

where the jump matrices J^ are also given in Figure |5] 
Denote u := e^'^*/^. As 2; — )■ oo with ± Imz > 0, we have 



^(z) 




1 lL±(/ + 0(^-i))5±e«(^^-), 



where L±, B± and 9(z; r) are defined by 

1 

B, 



L_ := 





B_ :-- 





1 








e"^ 

e-V 








eV 



(2.7) 



(2.8) 
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Q(^.^.^)(^^^Si0i{z;T),93iz;T),92{z;T)) for Im^ > 0, 
diag{92{z;T),63{z;T),6i{z;T)) for Im^; < 0, 



and the 9k are defined by 



for k = 1,2,3. 



(2.10) 



The expansion (2.7) for \l/(z) as z — )■ oo is valid uniformly for r in a 
bounded set. 



Denote by si, S2 the sectors 



si := <^^ G C 
S2 := < ^ e C 



ar£f2; G 



Stt 7r\ /tt Stt 






Around we then have the following estimate: 

Ui{z) 62(2;) e2(2;)\ 

^[z) = 0\tr{z) t2{z) t2[z)\ as 2^0, 

\ei(z) £2(2;) £2(2;)/ 



(2.11) 



with {tx{z),t2{z)) = < 



for /3 < 0, 
(l,log2;) for /3 = 0,2 e si, 

(log 2;, log 2;) for /3 = 0,2; G S2, (2.12) 
(2^,1) for /3>0,zGsi, 

Jl,l) for/3>0,zGs2. 



Note that the parameter /3 appears in the jump condition and in the 
behavior around 0. The dependence on r is only in the asymptotic condition 
as z — )■ 00. 

The RH problem for \E' has a unique solution. It is constructed out of 
solutions of the third order linear differential equation 



zq"{z) - l3q"{z) - Tq{z) + q{z) = 0. 



(2.13) 



This differential equation (2.13) has solutions in the form of contour integrals 

q{z)= f r^-''e^-^+'Ut, (2.14) 
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Ts 




Ts 





< — 





Figure 6: The contours Fi, F2 and F3 in the t-plane. The dashed hne denotes 
the cut of t-^-^ 



where F is an appropriate contour so that the integrand vanishes at the 
endpoints of the contour F. Define three contours Fi, F2 and F3 as in Figure 
|6l and define for z with Re z > 



?jU. 



-L 



t-^-^e^-^+^'dt, 



J = 1,2,3, 



(2.15) 



where we choose the branch of t ^ '^ with a cut on the positive real axis, i.e., 



The integrals (2.15) only converge for z with Rez > 0, but the functions 



Qj can be continued analytically using contour deformations. Branch points 
for the g-functions are and 00 and we take the analytic continuation to 
C \ (—00, 0], thus with a branch cut on the negative real axis. 

In the upper half-plane the unique solution for ^ is then given by, see 
for more details, 



m :-- 



/e2/3. 



qi 



g2/3^i^/ 






< arg z < J 



\^^2p^iq'i e/3^ig^' g^/y 



/e2/3. 



'qi + g2 e^'^^gs g2 



,\ 



e^l^^'q[ + g2 e^'^^g^ g^ 



g2/3«^// 



V 



4' 



Stt 



q'i e^^'q's q'ij 



_ „2/37ri. 



■gi + g2 - e^'^^'^gs 



e^^^'q[ + g2 - e^^^'q'^ 



aPm, 



e^"'g3 g2\ 
e^"*g^ q'2 



j<aigz< 4 , 



f <arg2<7r. 



(2.16) 



\e^f^^'q'l + q'i~e^^^'q'; e^^^g^ q'^J 
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and in the lower half-plane by 



^ := 



g2 e^'^'q^ 






e^^'q'^ -e^^^'q'l^ 



g^ + e'^^^'q[ 
\q'^ + e'^^WI 
(e^^^'qi + q2 



ef'^'q'l 



_„2/37rJ // 



^'i7 



93 



e^^^^gi + q'2 + q'z e^^'q'^ 



e^^^gs 



IV 



g2/3.i^// ^ ^^/ 



€ 



g/37ri^W 



-e^^^'qi\ 

.g2/3«^/ 
„2/37rJ„// 



9l/ 



■f <arg2;<0, 



-f <arg2;< -f, 



-TT < are 2; < 



" 4 • 



(2.17) 



3 Steepest descent transformations 

In this section we apply the Deift/Zhou steepest descent method to the RH 
problem for Y . For details of the various transformations we refer to [TS] . 

The Deift/Zhou method consists of a number of invertible transforma- 
tions reducing the original matrix-valued function y to a function R that is 
uniformly close to the identity matrix. Of main importance is the effect of 
these transformations on the expression for the MOP kernel i^„j^„2(x,?/), see 



(2.6), with diagonal multi-indices: 



Til 



n 



^2- 



(3.1) 



The first transformation Y ^ T normalizes the RH problem for Y at 
infinity: the new matrix-valued function T tends to the identity matrix at 
infinity. The transformation uses functions gi and g2 derived from a modi- 
fied equilibrium problem with logarithmic potentials. These ^f-functions are 
analytic on C\ (—00, 1]. See [TH Section 4.1] for precise details on the trans- 



formation Y \^ T. For the MOP kernel we obtain from (2.6) and formula 
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(4.4) in 

(g2n(gi,+{x)+92,+ (3::))\ 
. (3.2) 

The constants /i and I2 appear in the Euler-Lagrange variational conditions 
related to the equilibrium problem. 

The second transformation T 1— ;■ S* is the opening of the lenses, see |13t 
Section 4.2]. In this transformation the rapidly oscillating jump matrices 
along (a, 0) and (0, 1) that appeared after the transformation Y ^-^ T are 
turned into constant jump matrices, with the side effect of creating expo- 
nentially small jump matrices on the lips of the lenses. The expression for 
this transformation depends on the location of z. From now on we restrict 
ourselves to positive x and y, since in the cases x < or y < the expres- 
sions for Kn^x, y) are different but similar. Remark that we have Wiiy) = 



for y > 0. By (3.2) and formulas (4.14) and (4.5) in [13], we then get for 



X, y G (0, 1), (note that Wi{y) = 0): 



(_g-2n(gi,+(j,)+(;2,+ fe))^ T 

W2(y)e2"(s2,+(j/)+^2) 

(^2n{gi^+(x)+g2,+{x)) 
I . (3.3) 

iy2(a;)~"'"e~^"'^^2,+(a^)+'2) 

The next two transformations S ^ Rq and Rq ^ R are applied after 
constructing local and global parametrices. These parametrices are built as 
approximations to the RH problem for S, and one can think of Rq as the 
approximation error. The approximation around is not good enough for 
the further analysis, whence the need for a very last transformation Rq i— )■ R. 
We refer to [13| Section 4.4] for precise details and formulas. 

The local parametrix Pq around is constructed on a disk Uq centered 
in zero with shrinking radius ra"^/^, see formulas (4.68) and (4.70) of [13j. It 



involves the Angelesco model parametrix \E' described in Section 2^ Assume 
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that X and y are positive and belong to the disk Uq. Then S = RqPq and 
iserting thi 



inserting this in (|3.3|) we obtain after some calculations 

(-1 l)^+(n5/(y);r25r(y)) 



27ri{x — y) 



y"- 



X Eniy)-'Ro{y)-^Ro{x)E^{x)-^+ (nlf{xy,nh{x)] | | . (3.4) 



For X < or y < the formula (3.4) also holds, but with different row 



and column vectors, as given in the following table: 



y<0 



y>0 



X <0 




X > 



(-1 1 0) , 



(3.5) 



(-1 1), 



The functions / and r appearing in (3.4) were defined in 



Since we do 

not need to know their full expressions we will restrict ourselves to recalling 
their most important properties. 

These functions are both analytic in a neighborhood of 0, with /(O) = 0. 
They depend on a, and as a — )■ — 1 one has: 

/(O; a) = V2 + C(a + 1), r(0; a) = ^(a + 1) + 0{a + l)^ 



V2' 

see p^ formula (4.76)]. From this it easily follows that: 
Lemma 3.1. LetxER\ {0} and r eR. Put 



(3.6) 



X 



Xr. 



^nV 



-l + #. 



n2 



Then as n ^ oo we have 



n^f{xn;an) = x + O (n M 

n2T{Xn', an) = T + O iu^^ 

The O-term is uniform for x in a bounded set. 



(3.7) 

(3.8) 
(3.9) 
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Proof. This follows easily from (3.6) and the fact that f{z;a) and T{z,a) 



converge uniformly as a — > — 1 to f{z; —1) and t(z; —1] 

Also we have the following technical estimate: 
Proposition 3.2. Let x,y eM.\ {0}, r eM. and put 



n 



X 



Xr, 



V2nr 



Vn 



y 



V2- 



T 



V2nr 



1 
712 



Then as n ^ oo we have 

EniVn, anY^ RaiVnY^ Ro{Xn)En{Xn] a„) = / + O ( Ti^B 

The O-term is uniform for x and y in bounded sets. 



(3.10) 



(3.11) 



Assuming this proposition we can finish the proof of Theorem 1.1 



as 



follows. Let x,y > and let Xn,yn,0'n be as in (3.10). Then for n large 



enough we have that Xn,yn G Uq. Using (3.8), (3.9) and (3.11) in (3.4) we 
then obtain: 



1 



^/2nl 



-'^n.nyXjij i/ji) O'n) 



2TTi{x — y) 



(-1 l)^iy;T)-'^{x;T) io\+o(K 



(3.12) 



This proves Theorem |1.1| for x,y > 0. The right-hand side of (3.12) also 
provides us with an expression for ^^"^^{xyy-jT) in terms of the Angelesco 



parametrix \1/. For x < or y < the proof is completely similar and (3.12) 



holds, but with different row and column vectors, as given in (3.5) 



4 Proof of Proposition 3.2 



Proof. The proof of Proposition |3.2| is rather technical and requires deeper 
knowledge of the structure of the matrix functions Rq and En- 
From [131 Section 5.2] we have 



Ro{z)=R{z){I-Vn{z-a)), 



zeUn 



(4.1) 



where, see formula (5.13) of 

R{z) = 1 + 0(n-i/^) as n ^ oo, (4.2) 



uniformly in z. The function Ki in (4.1) is given by 



V^iz;a) = Z^iz;a)-^^^^ (4.3) 

A/ 

where Z„ and Z„ (a) play a role in the analysis of [13] . The function Z„ is 
analytic on f/o\{0} with a simple pole in 0, whose residue is Zn (a). Hence 
Vn is analytic on f/p. Additionally, from [131 equation (4.84)] it follows that 
Vn is of the form 

for a certain constant diagonal matrix D^. As a consequence we have the 
following basic property of V^: 

Vn{x\ a)Vn{y\ a) = 0, for every x,y e Uq, (4.5) 

Additionally, Vn{z; an) is uniformly bounded on the disk Uq. Since Uq has 
ius n~^^'^ t 
and therefore 



radius n ^1'^ this easily implies that l^(2;„;a„) = (9(n^/^) if z„ = (9(r2 ^/^ 



K(x„; a„) - K(yn; a„) = C(n ^). (4.6) 

From (4.1) and the fact that x„ = C(n~^/^), ?/„ = C(n~^/^), it follows that 

R{xin)-^R{xn) = I^o{n~'^. (4.7) 



Combining (4.7), (4.2), (4.5), and (4.6) we obtain 



i?o(l/n)-'i?o(x„) = / + 0[n-^). (4.8) 

The matrix-valued function En defined in formula (4.80) of [13j is analytic 
and invertible on a fixed neighborhood of 0. It depends on both a and n. It 
takes the form 



/3 riT(2;a) 



'n-2 
En[z\a) = n-^^e-"""^^E{z\a)\ 1 |, (4.9) 

ni 
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where E is independent of n, and E converges uniformly as a — ;■ — 1. Then 
E(y„; a„)-iE(x„; a^) = I + O (n~l) , (4.10) 



and combining with (4.9) yields 



En{yn;an) ^-E„(x„; a„) = / + O f 72 2 



(4.11) 



It also follows from (4.9) that 



n^2E^{xn; an) = O (n^j , n 2E„(i/„; a„) ^ = (n^j . (4.12) 



A straightforward combination of (4.8), (4.11) and (4.12), would yield an 



C(l) bound for the left-hand side of (3.11). We need to go a little deeper 
into the structure of E^ in order to improve the bound. 

From the computations in Section 6 of [13] it follows that En{xn', an) has 
the following behavior as ra — )■ oo: 

En{xn,an) = Cfi"-^ D^' I 1 (O l) (/ + O (^-i) ) , (4.13) 

where Ci is some non-zero constant, see [131 equation (6.12)]. From similar 
computations one finds for some non-zero C2: 



En{yn,an)'' = C2n^ io\{V2z 1 -1) D^ (^I + O (n-"^^y (4.14) 



From these expansions it follows that the leading terms of En{yn] an) and 
En{xn', an) cauccl out with Vn, see (4.4). In particular it holds by ( 4.4[ ), (4.6), 
( |ri3| and dili) ) that 



EniVn] an) ^[K(z/n; «„) - Vnixn] a„)]E„(x„) = 0{n ^) as u -^ 00. (4.15) 



Taking into account the extra estimate (4.15) we can perform the following 



computation (we do not write a„, and we use (4.1), (4.7), (4.5), (4.11), and 



20 



the fact that V^(x„) and Ki(|/„) are uniformly bounded) 






E^{y^,)-\I + Vn{yn)) 



l + 0[n—^ 



(/ - Vn{Xn))En{Xn) 



EniVn) 



I + Vn{yn)-Vn{Xn) + 0[n—e) E,(x„) 



which completes the proof of Proposition 3.2 



n 



5 Proof of Propositions 1.2 and 1.3 



In the previous section we derived expressions for K^'^s j^ terms of \Ef, see 



right-hand side of (3.12 ). In the present section these expressions will be used 



to obtain the explicit formulas stated in Proposition |1.2| and Proposition L3 
We follow here the computations in [HI Section 10.2], see also |23l Section 
8.2]. 

As stated in Section 12.21 the entries of \E' are solutions and derivatives of 
solutions to the third order differential equation 

zq"'iz) - (3q"iz) - rq'iz) + qiz) = 0. (5.1) 

The solution functions qj{z), j = 1, 2, 3 have explicit expressions 



Qji 



-13-3 



e* 2t- 



r+zt 



dt, 



(5.2) 



with contours F^ as shown in Figure [6j 

A first step towards the derivation of the double integral formula consists 
of finding similar explicit expressions for the entries of the inverse of \E'. This 



involves an adjoint differential equation to the one for the g-functions (5.1), 
and an associated bilinear concomitant, see e.g. [3], fT8] . 

The solutions of this adjoint differential equation have contour integral 
representations, and it takes some work to identify the appropriate contours. 
When the explicit expressions for the entries of \E'~^ are known we will be in 
a position to prove the double integral formula for ]K^°§. 
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The adjoint differential equation reads 

zr"'{z) + {P + 3)r"{z) - Tr'{z) - r{z) 



0. 



(5.3) 



The associated bihnear concomitant i3 is a differential operator defined on 
general twice differentiable functions q and r: 

B[q, r]{x, y) := yr{y)q"{x) - [(/3 + l)r(y) + yr'{y)\(({x) 

+ \vr\y) + (/3 + 1)r\y) - rr{y)\q[x). (5.4) 

In order to ease the notation we introduce 3 operators i3o, B\ and E^ 

Bor{z) := zr"{z) + {P + 2)r'{z)-Tr{z), 

Birlz) := -zr'{z)-{l3 + l)r{z), (5.5) 

B2r{z) := zr{z). 

So we can write the bilinear concomitant as 

B[q,r]{x,y) = B2r{y)q{x) + Bir{y)q{x) + Bor{y)q{x). (5.6) 

The central property of this bilinear concomitant is that if q solves the 



differential equation for q (5.1) and r solves the differential equation for r 



(5.3), then the bilinear concomitant evaluated on the diagonal B[q, r](a;, x) is 
a constant. This can be checked by taking the derivative of B[q, r]{x, x) with 
respect to x and plugging in the differential identities for q and r. 



The differential equation for r (5.1) allows solutions in the form of a 



contour integral: for an appropriate contour T the function 

as 



r(z] 



/e-^+^ 



(5.7) 



solves the differential equation (5.3). By an appropriate contour is meant a 



contour such that the integrand vanishes exponentially fast at the endpoints. 
In this case this implies that F must end and start in or at oo. 
Now we construct solutions rk{z), fc = 1, 2, 3 to (5.3) such that 



B[qj,rk]{z,z) 



Sjk j, k 



1,2,3, 



(5.8) 



where we put the cuts of the functions r^, if any, on the negative real axis. 
Using these relations we can find expressions for \l/~^. In each sector of the 
complex plane "^{z) is of the form 



qiiz) q2{z) qaizY 
^(z) = I q[{z) q'.iz) q',iz) 
.q'liz) q'^iz) q'iiz)^ 



■ c 



(5.9) 
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for some invertible matrix C. The inverse of \1/ is then given by 

fBoniz) B,n{z) B^niz) 

^{z)-' =C '■{ Bor^iz) B^z) B^r^iz) | . (5.10) 

\Bor3{z) Birsiz) B2r^{z) 

Indeed, the product ^ {z)^^^ {z) then equals 

(B[qi,ri]{z,z) B[q2,ri]{z,z) B[q^,ri]{z,z) 
C-^-\B[q,,r2](z,z) B[q2,r,]{z,z) B[qs,r2]{z,z) \ ■ C, (5.11) 

\B[qi,r3]{z,z) B[q2,r3]{z,z) B[q3,r3]{z,z) 



which by (5.8) is the identity matrix. 



The functions rk{z) can be chosen to be of the form 

rk{z) = ^ [ /e-5+5^-^^ds, (5.12) 

with the contours F^ determined by the following proposition: 

Proposition 5.1. Let Q G C be a non-empty open set. Assume that the 
following two contour integrals converge for z & Q: 

q{z)= I t-^'^'e^-^^^'dt, r{z)= [s^e'^^+^-''ds. (5.13) 

Also assume that the cuts of t^^^^ and s^ coincide, such that t^^^^s^ = t^^ 
for t = s. 

If r and r do not intersect, or only intersect in 0, then B[q,r]{z,z) = 
for z & Q. If they however intersect once transversally at a non-zero point, 
and if T meets T on the negative side of T, then 

B[q,r]{z,z) = 27ri, zeVL. (5.14) 

Proof. In order to simplify the notations we introduce the two phase functions 

d,{t):=''--^ + zt-l3\ogt, (5.15) 

er{s):=-- + —^-zs + {(5 + ?>)\ogs. (5.16) 
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6q and Or'. 



Using the definition of tlie bilinear concomitant (5.4) and tlie definitions 



of q{z) and r{z) (5.13) we can write tlie concomitant B[q,r]{z, z) of q and r 
as a double integral: 



z{s'^ + St + e) -{13 + 2)s - (/3 + l)t - r 



Mt)^e.{s)^g^l_ (5.17) 



'ter J s& 
The fraction in the above integral can be split up as follows: 

;z(s2 + st + t^) - (/3 + 2)s - (/3 + l)t - r 



deJs) 



1 ddq{t) s'^ + st + e 



(5.18) 



t3(t - s) ds s^{t - s) dt sH^{t - s) ■ 

Suppose that the contours F and F do not intersect, or only intersect in 



0. Then the three terms on the right-hand side of (5.18) remain bounded 
and we can write: 



B[q,r]{z,z) 



1 deJs) 



t&Js&t^it-s) ds 



'"^'^e'^^'Usdt 



+ 



1 



de,{t) 



Mt)Qer(s) ^^^^ 



terJsefs''^{t-s) dt 



s^ + st + t^ 



ter iser -^ '^ i't ~ s) 



(5.19) 



We can apply integration by parts to the first two integrals in (5.19). The 
boundary terms vanish. What remains is: 



B[q,r]{z,z) 



1 



+ 



0. 



ter Jser ^ [^ ~ -V 
1 



ter J ser 



S3(t-s)2 
t^- S^ 



ter Jser s t [t — s 






(5.20) 



Now suppose that F and F intersect in a point c 7^ 0, with F going from 
the right side of F to the left side. Then the contours can be deformed such 
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c-6 



c + i6 



c + 6^ 



c — iS 

r 



\ 



Figure 7: The deformed contours F containing the interval [c — i6,c + iS] and 
F containing [c — 6,c + 6]. 



that F contains the interval [c — i6,c + i6] while F contains [c — 6,c + S] for 
some 6 > 0, see Figure ul Define a new contour F^ by F^ := F \ [c — £, c + e] 



with < e < S. Then we have using (5.17): 



B[q, r]{z,z) = lim 



e-s>0+ 



ter Jser 



^(s^ + st + t^) -r 
+ (/3 + 2)s - (/3 + l)t 



,9,(t)^0,(s) 



S3t3 



(isfit. 



Since F^ and F do not intersect we can use the expression in (5.19) with F 



replaced by F"^ and apply integration by parts. However, now the boundary 
term in s remains: 



B[q,r]{z,z) 



lim , 



oSr{s) 



t 



'W dt. 



(5.22) 



=C+£ 



By continuity we have e^''*^'^^'^) = e^''-'^-' + 0{e) as e — > 0. Hence we obtain 

1 1 



%,r](z,^) = lim(e^'(^) + 0(£))y"i 



t — c + e t — c — e 



^Wrft. 



(5.23) 

Finally we deform the segment [c — i6,c + i6] C F into the union of a 
small circle around c + e and a semicircle with radius 6 around c, see Figure 
|8] The small circle around c + e will give us a residue contribution, while the 
contribution of the semicircle and the remaining part of F tends to zero as 
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c + iS {- 



•c + 5 



c — 15 • 



Figure 8: The dashed part of the contour F is deformed into a semicircle of 
radius 5 around c. Additionally a residue is picked up in c + e. The contour 
F is not shown here. 



e — >■ 0"''. Remark that the small circle around c + e has clockwise orientation. 
Then we have as e — )■ O"*": 



1 



t — c + e t — c — e 



3^'W dt = 27ri--^e^'('=+^) + 0(e). (5.24) 



Taking the limit e -^ 0^, and plugging the result into (5.23) we get 
B[q,r]{z,z) = 27ri^e^'(")e^'-(") = 27n, 



(5.25) 



where we used (5.15) and (5.16) 



n 

Choose contours F^ as in Figure |9J and let the functions r^ be defined by 

1 



rk[z) ■■-- 



2m ./r. 



s/^e"?+^-^^ 



ds, k = 1,2, 3, 



(5.26) 



where s^ has a cut on the positive real axis and s\_ = \s\^. Then it follows 



from Proposition 5.1 that the concomitant conditions (5.8) are satisfied. The 
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Figure 9: The contours Tk, k = 1,2, 3. The dashed hue denotes the cut of s^. 



integrals (5.26) define analytic functions on C with possible branch points in 



and oo and we put the branch cut on the negative real axis. 



Proof of Proposition L2. Having found explicit expressions for the entries of 



\1' ^ we can return to the expression for K'^"^ given in (3.12). The 1x3 and 



3x1 matrices pick out certain linear combinations of the Qj and r^, and we 
obtain 



with Qq and Tq given by 



[X 



y) 



B[qo,rQ]{x,y), 



(5.27) 



%{x) 



e^^^'qiiyx) 



for X > 0, 



27ri ' I e^'^'qi,+ {x) + e-^''*g2,+ (a;) for x < 0, 
,-2^TTi^_^^y^ for y > 0, 
g M-Ygd/) for y < 0. 

The functions gi and g2 have branch cuts on the negative real axis. Hence 



roiy) ■-- 



r2iy) 



e 



(5.28) 
(5.29) 



we have to specify that we need the positive boundary value in (5.28) in the 



case X < 0. Using contour deformation and a simple substitution in the 



integral one can see that go and tq are given by the contour integrals (1.15) 



and (1.16). 



n 



Finally we derive from (5.27) the double integral formula. 



Proof of Proposition 1.3 One can check that the derivatives of go and tq 



have the following expressions: 






2m' 



-x'-^\xf 



t-P+n-3^^-f-f^+t^^^ 



(5.30) 



ro 



27 



Sn) 



A (y) 



.l)n 



27li 



y'-\y\-'-' 



/+"e-^+& 



ds. 



(5.31) 



ro 



for n G Z"*" and x,y E M.\ {0}. The contours Fq and Fq are given in Figure 



3| From (|5.27|) and the definition of the bihnear concomitant (5.4) it follows 

1 



that we can write 



An'^iy - x)\y\ 



.P 



iero ^sero 



i/3+3 



yH"^ - (/3 + l)xyH + xy'^st 
+x'^ys'^ — (/3 + 2)x'^ys — rx'^y'^ 



TX X _|_J 






Introducing the phase functions rjg and rjr 



rig{t) :=-(3\ogt + 



TX 



X 

2t2 



r/,(.):=(/3 + 3)log.-^ + ;f^ 



dsdt. 
(5.32) 

+ 1, (5.33) 

(5.34) 



we can rewrite the double integral in (5.32) as 

[ J_ \ y^^"^ - (/3 + l)xyH + xy'^st 
tero Xefo ^ [ +^V - (/? + 2)x2ys - rx^y^ 



.n,{t)^vr{s) ^g^^_ (5 35) 



We have the following straightforward identity for the expression inside 
the square brackets: 

yh'^ - (/3 + l)xyH + xy'^st + x'^ys'^ - (/? + 2)x^ys - rx^y^ 
x^ys'^ drjr{s) xyH^ drjgit) 



+ 



xs — yt ds xs — yt dt 

4 S 4j.S 

X ys — xy r 3 3 X — y 



(5.36) 



(xs — yty ' " xs — yt 
The terms involving derivatives of the phase functions rjq and r]^ can be 
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simplified using integration by parts: 



tero Js&Fo s t xs — yt as 
1 xy'^t^ 



(5.37) 



ro J To 



sH^ {xs — yty 



^V,{t)^Vr{s) ^^^^_ 



(5.38) 



By plugging (5.37) and (5.38) together with (5.36) into (5.35) we find 



1 



sH^ 



tero JseVo 
= {x-y) 



yh^ - (/3 + l)xyH + xy'^st 
+x'^ys^ — (/3 + 2)x'^ys — rx'^y'^ 



oV,{t)^Vr{s) ^^^^ 






(5.39) 



From (5.32) we then conclude that K^'^s(x,|/; r) is given by 



K^-^(x,y;r) 



sgn(y) /■ f s^ 1 e * 2?f+* 
(2vr^)2 Zero Xefo Fxs - yt ^^-i^+s 



rfsrft. (5.40) 



n 
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